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Abstract. We have studied statistical properties of the values of the Wigner function 
W(x) of ID quantum maps on compact 2D phase space of finite area V. For 
this purpose we have defined a Wigner function probability distribution P(w) — 
(l/V) J 5(w — W(x))dx, which has, by definition, fixed first and second moment. In 
particular, we concentrate on relaxation of time evolving quantum state in terms of 
W(x), starting from a coherent state. We have shown that for a classically chaotic 
quantum counterpart the distribution P(w) in the semi-classical limit becomes a 
Gaussian distribution that is fully determined by the first two moments. Numerical 
simulations have been performed for the quantum sawtooth map and the quantized 
kicked top. In a quantum system with Hilbert space dimension iV(~ 1 /h) the transition 
of P(w) to a Gaussian distribution was observed at times t oc log.ZV. In addition, it 
has been shown that the statistics of Wigner functions of propagator eigenstates is 
Gaussian as well in the classically fully chaotic regime. We have also studied the 
structure of the nodal cells of the Wigner function, in particular the distribution of 
intersection points between the zero manifold and arbitrary straight lines. 
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1. Introduction 

The Wigner function (WF) pQ is an essential concept of phase space representation of 
quantum mechanics, namely it is a useful and faithful representation of a pure or mixed 
quantum state in terms of functions of canonical classical phase space variables. It 
has many applications in various branches of physics, in particular in quantum optics. 
However, WF cannot be interpreted as quantum phase space distribution as it can 
develop negative values, in particular due to well known oscillatory interference fringes 
following e.g. coherent wave-packet superpositions. 

WF also played an important role in the realm of quantum chaology [2J 
In particular, it has been conjectured that WF of stationary eigenstates of 
bounded dynamical systems in quasi-classical regimes localizes onto classically invariant 



Wigner function statistics in classically chaotic systems 



2 



components of phase space. For example, for classically regular phase space regions, 
like e.g. KAM tori, WF is supposed to become a Dirac delta function on a KAM torus 
whereas for a classically chaotic component WF is supposed to condense uniformly there. 
However, it is well known that this asymptotic behavior is to be understood in a very 
weak limit sense, namely WF is becoming uniform only after being integrated (with 
a very smooth test function) over non-small (or classical) region of phase space. On 
the other hand, on a smaller (quantum) scale, e.g. of Planck cell size, the phase space 
structure of WF of typical, say random or ergodic states, is very much unknown. This 
question is important for the understanding of decoherence and quantum stability with 
respect to system's perturbations as discussed recently |H 03 EE E| • 

In this paper we address the question of the structure of typical WF from a 
statistical point of view. We define and analyze the statistical distribution of values 
of WF of a given quantum state. In particular, we are interested in the relaxation 
of this distribution with time, when we start from initial coherent state, and in the 
corresponding time scale. Here we limit ourselves to classically fully chaotic and discrete 
time systems, namely the chaotic quantum maps, where the full phase space is the only 
topologically transitive ergodic component. 

We show that for chaotic quantum maps, the limiting WF value distribution is 
a Gaussian. The average value is fixed by normalization, while the second moment 
(or the variance) diverges in the quasi-classical limit (h — > 0), so we have roughly a 
symmetric distribution of positive and negative values of WF for random states. In 
addition, we show that for chaotic systems, the relaxation to equilibrium in statistics of 
WF happens on a short log H (Ehrenfest) time scale. Furthermore we show for chaotic 
systems, that the statistics of WF of typical chaotic eigenstates is the same as WF 
statistics of a random state, which is consistent with the established Berry's conjecture 
[2*Uj . Finally, we present statistical analysis of the structure of nodal cells of WF of 
chaotic or random states and show that it can be described by simplified models based 
on random trigonometric functions. 

In section 2 we outline the essentials of general Wigner- Weyl formalism and discuss 
two special cases of compact phase space, namely of toroidal and spherical phase space, 
which are later used in numerical examples. In section 3 we discuss WF value statistics 
of a random (or quantum chaotic) state. In secton 4 we discuss WF statistics on two 
numerical examples of quantum chaos, and analyze in particular the relaxation of WF 
statistics in quantum time evolution starting from an initial coherent state. In section 5 
we go beyond simple value statistics, and study also the phase space structure statistics 
of WF, such as its spatial correlation function and distribution of diameters of its nodal 
cells, for example. In section 6 we summarize our most important results and conclude. 

2. Weyl- Wigner formalism and the Wigner function 

Let us first review some of the essential general principles of Weyl- Wigner (WW) phase- 
space representation of operators (see e.g. |H]). As we are aiming at compact (finite) 
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phase spaces of chaotic maps, our discussion of WW formalism has to be a little bit 
abstract such as to be able to incorporate different topologies that are discussed later. 
We begin by considering a phase space x an d a Hilbert space H, such that a typical 
function A(x),x G X-> corresponds under quantization to a linear operator A over TL. 
The inverse map A — > A(x) is called a Weyl symbol, and if A — \if))(if)\ then the 
corresponding symbol is called the Wigner function (WF). The Weyl symbol can be 
formally constructed using a self-adjoint kernel operator lj(x) with the property 

ti{u(x)&(y)} = 5(x -y), lu(x)^ = u(x), (1) 

namely 

A(x)=ti{u(x)A}, (2) 
which generates also the inverse map 

A = J oIVA(x)Cj(x). (3) 

For the systems with classically compact (finite) phase space \ the corresponding 
Hilbert space is finitely, say iV dimensional. As a consequence, WW map (J2J) can 
be uniquely inverted only if the phase space is restricted to a finite set of N 2 points 
x' — i x nk] n,k = 0,1 ... N — 1}. Kernel is in this case defined in discrete points only, 
uj n k = uj(x n k) with the property 

tr{o) n fc^ m i} = 5 nm 5ki, (4) 

so that discrete WW map reads 

A nk = tr{u(x nk )A}, and A = 2j A n u nk . (5) 

nk 

The WF W(x) is defined as a phase space representation of the density operator p 
multiplied by a certain suitable normalization constant C 

W(x) = Ctr{uj(x)p}. (6) 

The constant C can be set for the convenience of a particular application, for example, 
usually it is set by the normalization of probability, J dxW{x) = 1. However, in this 
paper we are interested in the fluctuation of the WF so we determine the constant C by 
fixing the standard deviation a 2 = W 2 — W =1, where (...) denotes the average over 
phase space. We note that the kernel Cj is far from being completely specified by the 
property (0). In addition we need certain correspondence principles. For example, for 
the usual (non-compact) symplectic geometry x = (q,p) in d+d dimensional phase space 
the kernel reads &(q,p) = (27rh)~ d / 2 J dv e w ' v ' h \q + v /2) (q — v/2\. However, for compact 
geometries the WW map may not be unique. In such cases one may decide for the most 
reasonable or the simplest choice. In any case, as effective value of Planck constant 
vanishes any consistent choice should yield equivalent results in the semiclassical limit. 
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2.1. Wigner function on 2-d torus 

Many simple phenomena in the theory of bounded Hamiltonian (symplectic) dynamical 
systems can be demonstrated on a 2-d torus x = T 2 = [0,2vr) x [0,2vr). However, 
quantization of a torus is not trivial from mathematical point of view. Hence there are 
several different proposals for a WW formalism on a torus ^H]- Here we are following the 
reference U P to convenient scaling factors. Let the Hilbert space TC have dimension 
N, and the canonical bases, namely the position basis \n),n = 0,...,N — 1, and 
the momentum basis \k), k = 0, ... ,N — 1, satisfying (n\k) = iV _1//2 exp(— 2nink/N). 
The quantum phase space in our formalism is a discrete mesh of N x N points, 
x n k = (2irn/N,2irk/N). The WF function on such discrete phase space in the limit 
of iV — > oo mimics the continuous Wigner function [T]. Due to technical reasons we 
assume that N is odd. The WW kernel in this formalism is defined as 



]= exp{-27Tin'k/N)5{2l -2n + n')\l) (l + n'\ 



^>nk ~ 

^ N n'l 

where all indices run from — (N — l)/2 to (N — l)/2, which will be assumed whenever 
we refer to the quantized torus. For a pure quantum state ip WF is given by 

/ N 3 

W4n,k) = C t tr{u nk p} (7) 



22exp[-27rin , k/N]5(2l-2n + n')(l + ri\il>)(i/j\l) 

V ~ n',l 

where we define a fat Dirac delta function 5(1) as 

?/tn 1 / • /, /,r\ 1 sin(7r//2) 

2.2. Wigner function on a sphere 

As a second important special case we consider quantum mechanics of a spin J variable 
whose classical phase space can be identified with a unit sphere x — S 2 described by 
spherical angles 6 G [0,7r),ip G [0, 2tt). The WW formalism for SU(2) geometry was 
developed in ^3] with the kernel 

^0) = Tk q Yk q (x), x = (6, <p), 

kq 

where Y^ q are standard spherical harmonics and T^ q are multipole operators defined by 
E E (-l) J " m v / 2TTT ( 3 k J ,)\Jm)(Jm'\. 

m=-J m'=-J \ i / 

Symbols (_ J m k q ^,) are standard Wigner 3j symbols ^Zj. Using the kernel one defines 
WF as [12] 

23 q=k I 2 7+1 

W(x) = C s tr {u(x)p} = C S J2J2 GkMx), C s = JAtt^j-. (8) 

k=Q q=-k 
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The information about the quantum system state - density operator p is hidden in the 
coefficients 

G kq = tr{pfl q }. 

We note that SU(2) WF is defined continuously everywhere on a sphere. However, it is a 
superposition of finitely many spherical functions with coefficients Gk q which correspond 
to a discrete WF over a finite discrete mesh of points on some other compact geometries, 
like e.g. torus. 



3. Wigner function value statistics of a random state 

Quantum states of classically chaotic systems are usually associated to the so called 
ergodic or random wave functions. These random states can be constructed using very 
simple principles, like Berry's random plane wave superposition for billiards or 
chaotic analytic functions in Bargman or Husimi representation of Hannay j^j, see also 
In this paper we want to analyze the Wigner function of a random state and its 
statistical properties. To best of our knowledge this has not yet been attempted before. 
Starting from a fixed WF W(x) we define its value distribution P(w) as 

P(w) = 1 f 5{w - W{x))dx. (9) 



We assumed that phase space x is compact and has finite volume V. The averages with 
respect to the probability density P(w) will be denoted by (•••). The first moment of 
P(w) is fixed by normalization of Wigner function (given by constant C). The second 
moment is also fixed by purity of the state which implies 

dxW 2 {x) = C J dxW(x). (10) 

Since in this paper we are interested in fluctuation of WF it is convenient to determine 
the scaling constant C © by setting the standard deviation to one, a 2 = W 2 — W =1. 
This implies for the average values (first moments), using eqs. ()7I8|) : 

_ 1 1 

torus: W = , sphere: W = —=. (11) 

It may be instructive to define a relative standard deviation in units of average WF 
k = a/W. Then the above statement (fTTJ) says that the relative standard deviation 
k — ► oo of WF distribution diverges in the semi-classical limit iV — > oo (or j —>■ oo). 

We should note that the result on divergence of relative standard deviation is general 
and independent of the structure of the state. So, if we have a very non-uniform state, 
e.g. Gaussian wave-packet, then the statement says merely that the standard deviation 
is very large compared to the average value since we have almost all density concentrated 
around a small region. However, if we have a random (ergodic) state, then the WF value 
distribution should be transitionally invariant in phase space. Therefore, the values of 
first two moments of global phase space WF value distribution should also determine 
the average and standard deviation of local distribution around each phase space point. 
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In order to determine the entire WF value distribution we define an appropriate random 
WF model in a specific generic geometry, namely on a 2-d torus, although the result 
which will be obtained can be argued to be universal and geometry independent. We 
start with a random state in N dimensional Hilbert space 

N N 

1=1 1=1 

where q are random generally complex coefficients. It turns out (see e.g. ^H]) that for 
large N, coefficient q may be considered as independent complex Gaussian variables 
with variance (|q| 2 ) = 1/N (here and later (. . .) denotes averaging over an ensemble of 
random states (fT2")l). 

Translational invariance of WF value distribution implies that the distribution over 
phase space can be replaced with the distribution over ensemble of random states ([12)1 . 
so we choose to study distribution of WF at the most easily computable phase space 
point (0, 0) 



J (Vo| 2 + E^ + / >^')- ( 13 ) 



In the leading order 8(1 + 1') can be approximated by Kronecker Si+y so expression (fT3j) 
is a sum of N independent terms, the first |co| 2 is strictly non-negative, while the other 
have vanishing mean and finite fluctuation. Hence, due to central limit theorem, for 
large N ^ 1, the distribution of W 00 becomes Gaussian 

1 



P{w) = ^exp l-l(w - W) 2 j , (14) 

with the first and second moment which can be computed directly by ensemble averaging 
giving identical results to previous phase averaging (JTTJ), and a 2 = 1. 



4. Dynamics and relaxation of Wigner function value statistics 

In this section we will consider dynamics, namely statistics of WF of time evolving 
pure states. In particular we shall focus on systems with ergodic, mixing, and fully 
chaotic classical dynamics such that in the course of classical dynamics any non- 
singular initial classical measure relaxes to a uniform (micro-canonical) measure. When 
turning to quantum mechanics we pose a simple problem. Let us start with 'the most 
classical' initial state, namely with the coherent state (e.g. Gaussian wave packet): 
First, how does WF value distribution of a time evolving state relax into the Gaussian 
distribution (|T4*j) which characterizes the final ergodic (random) state? Second, what is 
the characteristic time scale of this relaxation process for a typical chaotic system and 
how does it scale with the value of an effective Planck constant? 
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4-1. Toy models of quantum chaotic dynamics 

In order to carry out the plan outlined above we need to define generic time evolution on 
the Hilbert spaces of the quantized torus or the quantized sphere, respectively, which are 
chaotic in the classical limit. This is not difficult as we simply consider popular models 
which have been widely studied in the literature, namely the quantized sawtooth map 
on the torus [18j, and the quantized kicked top [T§] . 

The quantized sawtooth map on the torus (classical counterpart is defined on a 
torus (q,p) G [0,27r] x [0, 2irL\) defined on the finite Hilbert space of dimension iV with 
the evolution operator 



where K is a kicking strength, T = 2-kL/N is a period of forcing and integer L 
(usually set to 1) measures the vertical size of toroidal phase space. We have introduced 
formal position and momentum operators with integer eigenvalues, namely n and rh 
respectively, satisfying n\n) = n\n), rh\m) = m\m). We should note that the classical 
sawtooth map is ergodic and uniformly hyperbolic for Kq > 0. As for the initial state 
we always choose a coherent state (which can be expressed in terms of Jacobi theta 
functions for toroidal phase space ^BJ) centered somewhere on the torus. 

As for the other system defined (classically) on a spherical phase space, the quantum 
kicked top of spin J has a unitary map Uk acting on a (2 J + l)-dimensional Hilbert 
space 



Here, parameter a is an angle of rotation in between the kicks and 7 is the kicking 
strength. In the paper we will consider the parameter values 7 = 7r/2 and a = 3, or 
a = 10, corresponding to classically mixed, or fully chaotic, phase space, respectively. 
Again, we prepare the system initially in the SU(2) coherent state [T3] placed somewhere 
on the sphere. 

For the purpose of illustration we first show the cascade of snapshots of WF 
starting from initial coherent state undergoing quantum dynamics with completely 
chaotic classical time evolution. The results are shown in figure for the quantized 
sawtooth map while similar-looking cartoons were obtained for kicked top. It is clear 
from the figure that after very short time (of the order of few kicks) the Wigner function 
relaxes into a universally looking distribution with roughly symmetric distribution of 
positive and negative values. We conjecture that this asymptotic distribution should 
be a Gaussian (fTH) as derived for a random model. Furthermore, we conjecture that 
the time-scale on which this relaxation takes place is just of the order of the so-called 
Ehrenfest time [21], that is the time needed that the initially localized wave-packet 
spreads over the accessible phase space t ehr = log ft/ A where A is a classical Lyapunov 
exponent. 




(15) 
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Figure 1. The dynamics of WF in the quantized saw-tooth map for parameters 
K = 0.5, L = 1 and dimension N = 101. The system is initially in the coherent state 
centered at (q,p) = (2/3tt,1/3it). WF at successive integer time steps is shown from 
left to right and from top to bottom. The color code bar of the Wigner function values 
is shown at the bottom of the figure. Note the significant contribution of negative 
values for longer times. 

4-2. Wigner function value distribution 

Firstly, we want to check the relaxation of WF value statistics P(w) starting from 
the initial coherent state. Using numerical experiments with the sawtooth map and 
the kicked top we confirm our expectation and obtain fast relaxation into statistically 
significant Gaussian distribution (^H). The cascade of distributions P(w) for several 
consecutive kicks in the sawtooth map is shown in figure 121 It is perhaps interesting 
to note that the relaxation goes through an intermediate distribution which seems to 
have exponential tails (see a snapshot at 10th kick). A very similar result is obtained 
also for the kicked top model shown in figure Ek- For comparison we show in figure 
stationary WF value distributions obtained from time evolution starting from a random 
state. In this case we obtain, as expected, for the initial state, and for time evolving 
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(a) (b) 

Figure 3. Time evolution of the Wigner function value distribution of the 
quantized chaotic kicked top for a — 10, 7 = ir/2 and J = 50. In (a) we take 
initial coherent state whereas in (b) we take initial random state. Consecutive kicks 
are show with different line styles and the full curve is a theoretical Gaussian H141 ). 



4-3. Relaxation time scale 



Secondly, we want to quantitatively characterize the deviations from a Gaussian 
statistics and thus to measure the time scale of relaxation. To this end we define the 
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excess e of WF value distribution P(w) 

e = ( w - W)ya 4 - 3. (16) 

Note, that e = for a Gaussian ()14j) and the size of e should roughly measure deviation 
from a Gaussian. As we have discussed above, we expect the relaxation process to take 
place within the Ehrenfest time scale 

t r ~logiV/A, (17) 

where A is an effective (average) Lyapunov exponent and N is a dimensionality of Hilbert 
space. This is certainly a lower bound to a relaxation time scale, but as we will show by 
numerical experiments, it also gives the right scaling and order of magnitude of the true 
relaxation time. For this purpose we choose only the quantized sawtooth map where 
much larger N is accessible so that the logarithmic scaling can be checked. 



N=81 
N=243 
N=729 
N=2187 
N=6551 




0.6 0.8 
t/log(N) 



Figure 4. Scaling of relaxation to Gaussian statistics in the quantized saw-tooth map, 
Kq = 0.5, L = 1, measured with the excess e of the distribution. The initial state is a 
coherent wave packet. Different dimensions N are indicated in the figure. 



In figure 0] we plot e as a function of time of WF value distribution starting from a 
coherent initial states, for different values of N over several orders of magnitude. Indeed 
we observe very clean transition from e ~ 0.65 to e = at around t ~ log N. In order to 
check also the dependence on the Lyapunov exponent \(K ) we have in figure El plotted 
the relaxation of excess e for several different values of chaoticity parameter K$ and 
fixed N. We have found that the transition time scale is proportional to A~ x therefore 
supporting the formula (fTTj) . 

Since e is only one number which cannot characterize the overall distribution P(w) 
we have made the following additional check. We know that for a random state in the 
limit N —>■ oo the probability of having negative value of WF goes to 1/2. Therefore we 
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Figure 5. Excess e as a function of time for the quantized saw-tooth map at different 
kicking strength Ko (see figure) and at fixed dimension N = 2187 and at L = 1. 
Lyapunov exponent is given by the formula: X(K ) — log((2+ii'o+((2+-Ko) 2 — 4) 1 / 2 )/2) 

HHj. 



propose to study the percentage of phase space with negative valued WF 

P-= f P{w)dw (18) 

J —oo 

as a function of time. The results of numerical experiment for the quantized sawtooth 
map are shown in figure E3 We see that P_ converges very fast to 1/2 on a time scale 
(i c ) proportional to log(iV) but significantly (around a factor of two) smaller than the 
relaxation time t T . We note that negative values of WF mean serious deviation from 
classical Liouville density, so our result indicates that quantum classical correspondence 
breaks significantly before the WF value distribution becomes actually close to a 
Gaussian. However, both time scales have the same scaling with N and A. 



4-4- Value statistics of Wigner functions of stationary eigenstates 
We note that the time evolution of a Wigner function of initial state 

|-0) = y^a fc l-0fc), a k E C, ^2\a k \ 2 = l, (19) 

k k 

can be written as 

W(x, t)=Y, ex P - W M, (20) 

kl 

W kl {x) =Ctr{£(z)|^ fc )(^|}, 

in terms of Wigner function basis W k i if \ipk) are eigenfunctions of the propagator Z7* 
with eigenvalues exp(itQ k ) . Therefore, stationary (or time averaged) properties of a 
WF W(x,t) are determined by the diagonal functions W kk , which are the usual WF of 
the eigenstates of U. It is therefore clear that for classically chaotic and ergodic systems 
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Figure 6. Percentage of phase space with negative WF as a function of scaled time 
for the quantized saw-tooth map, Kq = 0.5, L = 1 at several different dimensions N 
as indicated in the figure. 



one should expect almost all Wkk{x) to have the property of a WF of a random state. 
In order to confirm this conjecture we have computed the e& of WF value distribution 
of a complete set of eigenf unctions {Wkk(x), k — 1, 2, . . .} and analyzed its distribution. 
Indeed we found, as shown in figure for the quantized kicked top, that in the classically 
chaotic case almost all eigenfunctions have e ~ 0. This has been further quantitatively 
compared to Random Matrix Theory by computing excess distribution for Wigner 
functions of eigenvectors of a Gaussian orthogonal (GOE) random matrix. Indeed 
very nice agreement was found, meaning that also the number of Wigner functions 
with larger excess ('accidental scars') is within the statistical fluctuation predicted by 
Random Matrix Theory. In the same figure (|IJ) we also analyze excess distribution for 
mixed and regular classical dynamics, where the excess distribution has a nontrivial 
shape. 



4-5. Auto- correlation of the Wigner function 

Statistical description of WF at individual points in phase space is meaningful if 
statistical correlations between Wigner function values at different phase space points 
are small. Therefore it may be interesting to define and study the auto-correlation 
function of Wigner function, defined as 

1 f 

C(5x) — — I W(x)W(x + 5x)dx (continuous phase space), 
C(5n,5m) = ^ W nm W n+Sn ,m+Sm (discrete phase space). 

nm 

The dynamics of the auto-correlation and the corresponding Wigner function starting 
from the coherent state placed at (q = 2.1, p = 1.2) of quantized saw-tooth are shown in 
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Figure 7. Cumulative distribution of the excess e of the quantum kicked top 
eigenfunctions in the case of classically (almost) regular (a = 0.1,0.5), mixed (a = 3), 
and fully chaotic (a = 10) dynamics. To compare the classically chaotic case with 
random matrix theory, the result for GOE random matrix eigenfunctions is given in 
the inset in the blown up scale. 

figure |H| At small times, when Wigner function has still a clear non-random structure, 
the auto-correlation varies strongly with the displacement vector 5x = (5m, Sn). From 
autocorrelation function one can clearly observe the directions of stable and unstable 
classical flow. At later times C(Sx) becomes isotropic, and apart from the delta spike 
at the origin 5x = 0, equal a constant, namely C(Sx) ~ 5(8x) + const.. Indeed, for a 
random state we shall proceed to determine C(8x) exactly. 




-40 -20 20 40 -40 -20 20 40 -40 -20 20 40 
kick:2 kick:5 kick:10 



Figure 8. The density plot of autocorrelation C(8x,5y) (below) and Wigner function 
W (x,y) (above) in time evolution of the sawtooth map at dimension TV = 101 and 
K = l,L=l. 
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We assume a random wave hypothesis, namely that a random (time-dependent) 
state can be written as \ip) = J2 n c n \n), where c n are Gaussian random independent 
complex coefficients, satisfying (c*c m ) = 5 nm , (c n c m ) = 0, where (. . .) represents the 
average over an ensemble of states, or equivalently, a time-average with a random initial 
state. The auto-correlation function can thus be expressed as 

C(Sn,Sm) = j^-j\a(5n,5m)\ 2 , (21) 

a(5n,5m) = exp(— \2itl5m/ N)cic* l+5n . 
i 

Averaging over Gaussian random variables be performed straightforwardly, 

yielding 

N 2 

(C(5n,5m)} = — - - 6 Sn>0 S Sm> o + W . 

This theoretical prediction is in good agreement with numerical experiment, as shown in 
figure El We also estimate (temporal) fluctuations of C(5x ^ 0) by standard deviation 
o 2 c = (C 2 ) — (C) 2 , which can be computed directly using Wick pair contractions in the 
variables c n 

o-c = (C) = j^-j = W 2 if (Sn, 5m) ^ 0, 
in the asymptotic regime of high dimension N ^> 1. 

2 
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Figure 9. Scaled local time average auto-correlation function C versus the radial 
distance r — y/ (Sx) 2 + (5y) 2 for different dimensions N of the quantized sawtooth 
map for Kq = 1 and L = 1. 



5. Wigner function phase space structure statistics 

In previous sections we have been investigating the value statistics of the Wigner 
function. WF possesses hills and valleys as compared to the reference "altitude" - 
offset. Here we discuss 2d compact phase space where these structures are supported by 
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2d nodal cells, that in general have very rich topology. The nodal cell sizes, especially 
the sub-Planck size structures, have recently been brought in the connection with 
decoherence |H E2] • We are now interested in the statistical properties of WF phase 
space structures such as the nodal cell size and the amplitude of oscillations (hills and 
valley) from a given offset value. For the offset we choose the phase space mean value 
of WF (offset = W) as this seems to be the most natural choice. However, since the 
relative mean value divided by the fluctuation goes to zero, W /aw — > in the limit 
N — > oo, we argue that the results are asymptotically just the same as if one considers 
nodal cells with respect to zero offset. 

In this work we would like to obtain some general results on WF nodal cell statistics 
of random (or chaotic) states. In our numerical analysis we consider time-evolving states, 
starting from initial coherent state, for times t > when the statistics is stationary 
and the structure of the state is the expected to be the same as that of a random state. 
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Figure 10. The sketch of ID phase space structures of the intersected WF. 

In order to simplify the picture we have considered the statistical properties of 
intersections between nodal cells (bounded by the curves W(x) — W = 0) and some 
arbitrary (random) straight lines in phase space, i.e. we consider one-dimensional 
projections of two-dimensional WF. Let the straight line R(t) G x be parametrized 
by a real variable t. In the case of spherical phase space, the set of all straight lines 
consists of all big circles, while in the case of toroidal phase space we consider for 
simplicity only closed straight lines, i.e. the two sets of irreducible circles specified by 
either fixing position or fixing momentum. 

The central object studied here is the so-called Wigner function on a line 
(abbreviated WFL) with an offset value subtracted and assuming that the circumference 
of the circle is equal to 2tt: 

W i>tR (t) = W i> [R(t)]-W, R(t)e x , te[0,24 (22) 

Since the statistical properties of a Wigner function of a random or ergodic quantum 
state are invariant under phase space translations, the parameters of the line (circle) 
R(t) can also be chosen as random. In other words one may average over lines R(t) 
in a uniform way, such that an ensemble of lines R(t) uniformly covers the phase 
space with respect to an ergodic invariant measure. We shall investigate the following 
statistical distributions defined with respect to W^^it): (i) distribution dV(s)/ds of 
spacings s = t n+ i — t n between adjacent zeros t n of W^r{P) which may also be called 
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'diameter distribution of Wigner nodal cells', (ii) distributions of amplitudes dP(A)/dA, 
i.e. local maxima (hills) and local minima (valleys) between each pair of adjacent zeros 
as illustrated in figure 1X0} and (iii) joint distribution of spacings s and the corresponding 
amplitudes A, dV(s, A)/dsdA. 



5.1. Random model of a Wigner function on a line 

We would like here to propose a simple statistical model which reproduces the properties 
of a random Wigner function on a line. We start by a simple ansatz expanding WFL 
into the Fourier modes 

M 

W (t) = u + u q cos(gt) + v q sin(gt) , (23) 

9=1 

where coefficients u q and v q are some Gaussian (asymptotically as N — > oo, statistically 
independent) random variables with zero mean (u q ) = 0, (v q ) = and prescribed 
variances (uV) = 0(iV _1 ), (v?) = 0(A r ~ 1 ). The effective number of Fourier modes 
M is usually of the same order than the dimension of the Hilbert space N. The ansatz 
(|23|) shall be proven separately for the case of WF on the torus and on the sphere, and 
also the expressions for the variances of coefficients u n and v n shall be computed. We 
will then use our statistical model to compare with exact numerical calculations in our 
two model systems. 



5.1.1. Random model on the torus We expect that statistical structure of random WF 
to be the same in both simplest sets of directions, namely for fixed position or fixed 
momentum. We limit the analysis which follows to the circles of fixed position. We 
start from a 'random' state written in the position basis {|ra)} n =i,...,7v of the Hilbert 
space of dimension N, as 

JV 

l^> = ^c», c n eC, (24) 

n=0 

with uncorrelated complex Gaussian coefficients (in the asymptotic regime N — > oo) 
specified by 

(cn) = 0, (c* m c n ) = — , (c m c n ) = 0. (25) 

The expression of WF (J2J) at some fixed position n can be rewritten in terms of 
continuous momentum variable t = 27im/N e [0, 27r) 

m 

W^ n {t) = Z i( n ) exp(igt), M = {N- l)/2. (26) 

q=-M 

Complex Fourier expansion coefficients Z q are expressed in terms of state coefficients c n 

N - 

V N - 1 l= _ M 
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Now, the Fourier modes in (J23j) are a simple bilinear functions of the coefficients c n , 
namely 



u = Z -W, 



Un 



Za + Z* 



i(Z n 



q = l,2...M. (2* 



By means of the central limit theorem one can argue in the asymptotic regime of large 
N, due to summation of statistically independent terms in eq. (|27jh that u n and v n 
become (independent) Gaussian variables of vanishing first moments (u n ) = 0, (v n ) = 0. 
Straightforward but tedious calculation gives for the second moments 

2 1 

{u q u q ,) = (v q v ql ) = — — -5 q>ql , (v q u q/ ) = 0, (u 2 ) = — — -. (29) 



Therefore, the variances of Fourier coefficients are asymptotically, as iV — > oo 
independent of the mode number q. This has been also verified by means of numerical 
simulation of WFL of random states, where c n have been generated using a suitable 
random number generator. The results are shown in figure ITTk. where one sees excellent 
agreement with the prediction ([29)1 . 
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Figure 11. Variances of the Fourier coefficients of WFL of a random state: (a) in 
case of toroidal phase space, and (b) in case of spherical phase space. 



5.1.2. Random model on the sphere All big circles on the sphere are statistically 
equivalent with respect to WF (jHJ) of a random state. Thus we choose to consider 
the equator 9 = n/2 for simplicity. Using explicit expression for the spherical harmonics 
Y kq 

Y kq {x) = N kq P q k {cos{6)) exp(z#), (30) 

where N kq are normalization constants and are generalized Legendre polynomials, 
we have explicit expression for the WF on the equator 

2J 

W4<j))= ^exp(z#), (31) 
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in terms of complex Fourier coefficients Z q 

Z q = C s 2j C m,C m tK q mTn ,, 
m,m'=—J 
2J 

Km* = E N kq P«(0)(Jm'\Tl q \ Jm) e R. 

k=\q\ 

The coefficients Z q are bilinear forms of wave coefficients c m with the Gramm matrix 
K q mm ,. The latter is difficult to evaluate explicitly. Using the symmetry of coefficients 
Z q = Z*_ q) we rewrite WFL as 

2J 

W^{4>) - W = u + ^2 u i cos (#) + v q sin(g0) , (32) 

where real Fourier coefficients are again simply related to complex coefficients Z q by 
([28)1 . The WF random model of spin J system has M = 2 J modes. In the semi-classical 
regime (large J) we can again show by means of central limit theorem that u q and v q 
should be Gaussian distributed. Taking into account the detailed properties of terms in 
Z q we can find the following statistical properties of u q , v q and u&. 

(Ug) = (Vg) = (U ) = 0, (UgV q/ ) = 0, 

(u q Uq') = (v q V q r) = (T q d q , q >, q = 1, . . . , 2 J. 

The standard deviations of coefficients o~ q are essential for further discussion, but they 
are quite difficult to compute analytically. We therefore obtain them by numerical study 
of WFL for a random state. The results of numerical simulation are shown on the figure 
fTTlb. For large J, the results seem to be perfectly fitted by the semicircle law 




< = iz\i l ~{^f) . ff=l,...,2J, (33) 
which is conjectured to be the correct semi-classical limit. 



5.2. Numerical study of the structure statistics of Wigner function 

Here we report on numerical simulations of the structure statistics of WFL of time- 
dependent states in two classically chaotic quantum systems, namely of the quantized 
sawtooth map at K = 10 on the torus, and quantized kicked top at a = 10,7 — 7r /2 
on the sphere. In addition, we simulate for comparison the corresponding statistics for 
the appropriate random model discussed in the previous paragraphs, since these seem 
to be impossible to be expressed analytically. 

Results of our numerical studies of structure statistics dP / ds, dV/dA and dV/dsdA, 
are shown on the figure Numerical experiments are done so that they correspond to 
WF random model with M = 50. We can see a good agreement between the random 
model and numerical measurements on real dynamical systems. The latter are performed 
by time averaging over the evolution of an arbitrary initial state. We should stress that 
the numerical results are identical if we instead consider Wigner function of a random 
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(ergodic) state. With this results we have again shown, indirectly, that the assumptions 
in the derivation of the random model are correct. In the case of a spherical phase 
space there is a small deviation of the spacing distribution dV/ds around the peak of 
distribution due to imperfect approximation of the mode variances with the asymptotic 
semicircle formula (|HBJ). 

The spacing distribution between neighboring zeros dV/ds of WFL (shown 
in figures 1121 ) gives some information on typical intersections of the nodal cells, 
whereas amplitude distribution dV/dA (shown in figures IT2*tI) gives information on 
the distribution of the heights in such intersected filaments. However, we note an 
interesting observation, namely that the statistics of the structure of positive nodal 
cells (W — W > 0) are identical to the statistics of negative nodal cells (W — W < 0) 
in the asymptotic regime N — > oo. This is consistent with an asymptotic symmetric 
Gaussian distribution of Wigner function values (J14)) . The spacing distribution exhibits 
periodically spaced peaks, with period ir/M, and exponential tail with non-universal 
exponent for large s 

As the physical domain of Wigner function on a torus is in fact a discrete 
mesh of N x N points, we also consider statistics of a proper discrete WFL, namely 
w(m) = W(27rm/N), and the discrete spacing distribution P s which in fact measure the 
distribution of clusters of points where w(m) has constant sign, and the corresponding 
amplitude distributions. These are shown in figure IT2*k . Interestingly, we find strong 
numerical evidence of the following exponential distribution 

P s = 2- s . (34) 

On the other hand, this is in fact just the distribution of the lengths of clusters of 
repeated uncorrelated binary events. This is consistent with the statement that the 
values of WF of a random state on a physical quantum mesh of N x N phase space 
point are indeed uncorrelated as suggested already by the auto-correlation function. We 
have also investigated the joint spacing-amplitude distributions dV(s, A)/dsdA, shown 
in figure fTBl Again the results inferred from the simulations of the random model match 
the simulations of WFL. We found that distributions dV/dsdA have a simple shell- 
shaped form, indicating strong correlations between the spacing and the amplitude of 
the WF intersections. In fact one can give an upper bound on the maximally allowed 
amplitude A at a given small spacing s, namely this can be estimated from the second 
derivative (curvature) of the random model ansatz (j2HJ). 

6. Summary and discussion 

In this paper we have proposed to study statistical properties of Wigner functions of 
random pure quantum states, either eigenstates or time-dependent states of classically 
chaotic systems. We concentrated on the properties of quantum maps with 2d classical 
phase space. We have shown, using simple general arguments, that the Wigner function 
value distribution of a random state should tend to a Gaussian, in the semi-classical 
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Figure 12. The average distribution of spacing s (I) and amplitude A (II) of WF 
along a random line in phase space: (a) discrete torus (see text for explanation), (b) 
continuous torus, (c) sphere. In the cases (b) and (c) the results of WFL simulations 
are compared with results obtained by simulations of the random model for the specific 
geometry. In all examples we consider the case with M = 50 Fourier modes. 



limit, which is centered around zero value, so the probabilities to have a negative or 
positive value become equal in the limit. In other words, the standard deviation divided 
by the mean value of the Wigner function diverges as a/W ~ y/N in the semi-classical 
limit N — > oo. In addition, we have analyzed the structure and phase-space correlations 
of Wigner functions of random states. In particular, we have shown that the auto- 
correlation of the Wigner function becomes a delta function in the semi-classical limit, 
and that the nodal cells have typical structures on (sub-Planckian) scale 5q, 5p ~ h. 
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Figure 13. The average joint distribution of amplitude A and spacing s between two 
successive zeros of WFL: (a) continuous torus, (b) sphere. For more details see caption 
of figure E| The density plot with grayness scales (inicated in the horizotal bar) refers 
to the simulation of true WFL, whereas isodensity contours refer to random models as 
discussed in the text. 



We believe that our results may shed some new light onto the related studies of 
decoherence [3J 122]- in particular, one now expect certain properties of the Wigner 
function of random states to be manifestly non- classical. For example, the fidelity of 
two initially equivalent states which undergo two slightly different time-evolutions has 
been found to behave non-classically for times larger than the — log h~ Ehrenfest time, 
i.e. when the states become effectively random, whereas classical behavior of fidelity 
has been recovered for shorter times [Sj. We stress that we have in this study only 
considered non-autonomous (time-dependent, e.g. kicked) quantum systems, where 
Wigner functions of ergodic states are expected to occupy the entire classical phase 
space. On the other hand, one may ask similar questions about the statistics of Wigner 
functions on, or close to, energy surfaces of autonomous Hamiltonian systems. This is 
a subject of a forthcoming publication [25]. 

Acknowledgments 

Useful discussions with A. Backer, G. Veble and M. Znidaric, as well as the financial 
support by the Ministry of Education, Science and Sport of Slovenia are gratefully 
acknowledged. 

References 

[1] Wigner E 1932 On the Quantum Correction of Thermodynamic Equilibrium Phys. Rev 40 749-59 
[2] Berry M V 1977 Semi-classical mechanics in the phase space: A study of Wigner's function 

Phil.Trans. R. Soc. A 287 237-71 
[3] Voros A 1976 Annls. Inst. Poincare A 24 31-90, Voros A 1977 Annls. Inst. Poincare A 26 343-403 
[4] Zurek W H 2001 Sub-Planck structure in phase space and its relevance for quantum decoherence 

Nature 412 712-5 



Wigner function statistics in classically chaotic systems 



22 



[5] Prosen T and Znidaric M 2002 Stability of quantum motion and correlation decay J. Phys. A 35 
1455 

[6] Karkuszewski Z P, Jarzynski C and Zurek W H 2002 Quantum Chaotic Environments, the Butterfly 

Effect, and Decoherence Phys. Rev. Lett. 89 170405 
[7] Jacquod Ph, Adagideli 1 and Beenakker C W J 2002 Decay of the Loschmidt Echo for Quantum 

States with Sub-Planck-Scale Structures Phys. Rev. Lett. 89 154103 
[8] de Groot S R and Suttorp L G in "Foundations of Electrodynamics" , Chapter VI (North-Holland, 

Amsterdam 1972) 

[9] Hannay J H 1998 The chaotic analytic function J. Phys. A 31 L755-L761 

[10] Nonnenmacher S and Voros A 1998 Chaotic Eigenfunctions in Phase Space J. Stat. Phys 92 431-518 
[11] Agam O and Brenner N 1995 Semiclassical Wigner functions for quantum maps on a torus J. 
Phys. A 28 1345-60 

[12] Dowling J P, Agarwal G S and Schleich W P 1994 Wigner distribution of a general angular- 
momentum state: Application to a collection two-level atoms Phys. Rev. A 49 4101-9 

[13] Arecchi F T, Courtens E, Gilmore R, Thomas H 1972 Atomic coherent states in quantum optics 
Phys. Rev. A 6 2211-37 

[14] Agarwal G S 1981 Relation between atomic coherent-state representation, state multipoles, and 

generalized phase-space distributions Phys. Rev. A 24 2889-96 
[15] Kasperkovitz p and Peev m 1994 Weyl- Wigner formalisms for toroidal geometries Ann. Phys. 230 

21-51, Luis A, Pefina J 1998 Discrete Wigner function for finite-dimensional systems J. Phys. 

A: Math. Gen 31 de Almeida O, Rivas A M F 1999 The Weyl representation on the torus Ann. 

Phys. 276 223-56 

[16] Leboeuf P and Voros A 1990 Chaos-revealing multiplicative representation of quantum eigenstates 

J. Phys. A: Math. Gen. 23 1765-74 
[17] Edmonds A R 1974 Angular momentum in quantum mechanics (Princeton, Princeton University 

Press) 

[18] Benenti G, Casati G 2001 Sensitivity of Quantum Motion For Classical Chaotic Systems, 

arXiv:quant-ph/0112060/l 
[19] Haake F 1991 Quantum signatures of chaos (Berlin: Springer) 

[20] Berry M V 1977 Regular and irregular wavefunctions J. Phys. A: Math. Gen. 10 2083-91 
[21] Bogomolny E, Bohigas O, Leboef P 1996 Quantum Chaotic Dynamics and Random Polynomials 
J. Stat. Phys. 85 639-79 

[22] Jordan A, Srednicki M 2001 Sub-Planck Structure, Decoherence, And Many-Body Environments, 

arXiv:quant-ph/0112139/l 
[23] Kuznetsov P L, Statonovich R L Tikhonov V I 1954 On the duration of the exceedences of the 

random function J. Tech. Phys. Moscow 24, 103-112 
[24] Berman G P and Zaslavsky G M 1978 Condition of stochasticity in quantum non- linear systems 

Physica 91A 450 
[25] Prosen T and Backer A 2002 in preparation 



